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Abstract — The dual of an entanglement-assisted quantum 
error-correcting (EAQEC) code is the code resulting from ex- 
changing the original code's information qubits for its ebits and 
vice versa. As an introduction to this notion, we show how 
entanglement-assisted repetition codes and entanglement-assisted 
accumulator codes are dual to each other, much like their classical 
counterparts, and we give an explicit, general quantum shift- 
register circuit that encodes both classes of codes. We later show 
that our constructions are optimal, and this result completes 
our understanding of these dual classes of codes. We obtain 
linear programming bounds for EAQEC codes, by exploiting 
general forms of these dualities and corresponding MacWilliams 
identities. We establish the Gilbert- Varshamov bound and the 
Plotkin bound for EAQEC codes, and all of these bounds allow 
us to formulate a table of upper and lower bounds on the 
minimum distance of any maximal-entanglement EAQEC code 
with length up to 15 channel qubits. Finally, we provide an upper 
bound on the block error probability when transmitting maximal- 
entanglement EAQEC codes through the depolarizing channel. 

Index Terms — quantum dual code, entanglement-assisted 
quantum error correction, MacWilliams identity, linear 
programming bound, entanglement-assisted repetition codes, 
entanglement-assisted accumulator codes 



I. Introduction 

The theory of quantum error correction is the theory under- 
pinning the practical realization of quantum computation and 
quantum communication flT] , 0, 0, 0, 0, 0. Quantum 
stabilizer codes are an extensively analyzed class of quantum 
error-correcting codes because their encoding, decoding, and 
recovery are straightforward to describe using their algebraic 
properties Q, 0, 0, ifTUl . In particular, a quantum code 
designer can produce quantum stabilizer codes from classical 
binary and quaternary self-orthogonal codes by means of the 
CSS and CRSS code constructions, respectively OH, QZ), 0. 
0. 

Entanglement-assisted quantum error correction is a 
paradigm for quantum error correction in which the sender 
and receiver share entanglement before quantum communi- 
cation begins [13]. An [[n, k,d;c\] EAQEC code encodes k 
information qubits into n channel qubits with the help of c 
pairs of maximally-entangled Bell states. The code can correct 
up to L^f^J errors acting on the n channel qubits, where d is 
the minimum distance of the code. Standard stabilizer codes 
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are a special case of EAQEC codes with c = 0, and we use 
the notation [[n, k, d]] for such codes. 

Bowen proposed the first EAQEC code fl4l . by replacing 
the two ancilla qubits of the three-qubit bit-flip code with 
two pairs of ebits. He then applied a unitary transformation 
on the code space to make the EAQEC code equivalent to 
the well-known five-qubit code [15]. Fattal et al. established 
a technique for handling entanglement in the stabilizer for- 
malism |[T6l . Brun, Devetak, and Hsieh then devised the 
entanglement-assisted stabilizer formalism and showed how 
to transform any [n, k, d] classical quaternary codd^ into an 
[[n, 2k — n + c, d; c]] EAQEC code, where c depends on the 
properties of the classical code [17|. Lai and Brun further 
explored the properties of EAQEC codes and proposed an 
optimization method to find optimal EAQEC codes that cannot 
be obtained by the aforementioned construction [ 1 8 1 . An 
[[n, k, d; c]] EAQEC code is optimal in the sense that d is the 
highest achievable minimum distance for given parameters n, 
k, and c. 

In this paper, we define the dual of an EAQEC code 
as the code obtained by exchanging the information qubits 
and the ebits of the original EAQEC code. We establish 
the MacWilliams identity for EAQEC codes, which we then 
exploit to construct a linear programming bound on the 
minimum distance of an EAQEC code. We can now show that 
several code parameters proposed in Ref. lTT8l are optimal, by 
exploiting the linear programming bound developed here. 

We establish the Gilbert- Varshamov bound for EAQEC 
codes, proving the existence of EAQEC codes with certain 
parameter values. We then apply the encoding optimization 
algorithm from Ref. |fl8l to find good EAQEC codes with 
maximal entanglement for n < 15. All of these EAQEC 
codes have minimum distance greater than or equal to that 
given by the Gilbert- Varshamov bound. Lai and Brun recently 
proposed a family of [[n, 1, n; n — 1]] entanglement-assisted 
repetition codes for n odd [18|; in this paper, we propose a 
family of [[n, 1, n — 1; n — 1]] entanglement-assisted repetition 
codes for n even, which we prove here to be optimal, and 
thus completing the family of entanglement-assisted repetition 
codes for any n. We also derive the quantum version of the 
Plotkin bound |[T9l . which is tight for codes with small k 
and maximal entanglement. Combining the results of linear 
programming bounds with the existence of EAQEC codes, 
we give a table of upper and lower bounds on the highest 
achievable minimum distance of any maximal-entanglement 

'An [n,k, d] classical linear code over a certain field encodes k information 
digits into n digits, where d is its minimum distance. 
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EAQEC codeQfor n < 15. 

The weight enumerator of a classical code gives an upper 
bound on the block error probability when transmitting coded 
bits through a binary symmetric channel 11241 . (25]. Since 
maximal-entanglement EAQEC codes have many similarities 
with classical codes |fl~8), 11231 . we can find an upper bound 
on the block error probability of transmitting coded quantum 
information through the depolarizing channel, and this deriva- 
tion is similar to the classical derivation 1241 . We also exploit 
this result to find an upper bound on the expected block error 
probability when decoding a random maximal-entanglement 
EAQEC code. 

We organize this paper as follows. We first review the 
notion of a dual code in classical coding theory and recall the 
MacWilliams identity for quantum codes. We then give the 
definition of the dual of an EAQEC code, and we explain this 
notion with the example of the dual repetition and accumulator 
EAQEC codes. We follow the terminology and notations of 
EAQEC codes used in [ 18]. For details, we point the reader 
to Refs. flTJ], AS). The MacWilliams identity for EAQEC 
codes and the linear programming bound for EAQEC codes 
are derived in Sections IV and [V] In Section VI we begin 
with the Gilbert- Varshamov bound for EAQEC codes. Then we 
describe the construction of [[n, 1, n — 1; n— 1]] entanglement- 
assisted repetition codes for n even and prove other results 
about the existence of EAQEC codes. We finish this section 
with a table of upper and lower bounds on the minimum 
distance of any EAQEC codes with maximal entanglement 



for n < 15. Section VII details an upper bound on the block 
error probability under maximum-likelihood decoding, and the 
final section concludes with a summary and future questions. 

II. Preliminaries 

In classical coding theory, a well-established notion is that 
of a dual code. Suppose that C is an [n, k] linear code over 
an arbitrary field GF(q) with a k x n generator matrix G and 
a corresponding (n — k) x n parity check matrix H such that 
HG T = 0( n -k)xk- The transpose of the above matrix equation 
holds, i.e., GH T = Ok X (n-k)i revealing that H can play the 
role of a generator matrix and G can play the role of a parity 
check matrix. The dual code of C is the [n, n — k] linear code 
C with H as a generator matrix and G as a parity check 
matrix. 

Suppose the weight enumerator of a linear code C is 



A(z) 



w=0 



An 



where A w is the number of codewords in C of weight w. The 
generalized MacWilliams identity gives a relation between the 

2 One might wonder why we are considering EAQEC codes that exploit the 
maximum amount of entanglement possible, given that noiseless entanglement 
could be expensive in practice. But there is good reason for doing so. The 
one-shot father protocol is a random entanglement-assisted quantum code 
1201 . 1211 . and it achieves the entanglement-assisted quantum capacity of a 
depolarizing channel (the entanglement-assisted hashing bound [221, |14'|) by 
exploiting maximal entanglement. Furthermore, there is numerical evidence 
that maximal-entanglement turbo codes come within a few dB of achieving 
the entanglement-assisted hashing bound 1231 . 



weight enumerator A(z) of C and the weight enumerator B(z) 
of its dual code C 1 - (H): 

B ^=U 1+{q - 1)z)nA{ ^^ ] - (1) 

Consequently, we can determine the minimum distance of the 
dual code C from the weight enumerator of C. 

A good example of this duality occurs with a repetition code 
and an accumulator code. The [n, l,n] binary repetition code 
C has the lxn generator matrix 



G = [1 1 1 • • • 1 1] 
and an (n — 1) x n parity check matrix 



(2) 



H = 



1 1 
1 1 









1 1 



Its weight enumerator is A(z) = 1 + z n . Flipping the roles of 
these matrices gives an [n,n — 1,2] binary accumulator code 
with the following weight enumerator: 

B(z)=h(l + z)» + (l-zr) 



E 



Prior to this paper, the notion of a dual code did not exist 
in quantum coding theory. However, the quantum analogue of 
the MacWilliams identity does exist ll26l . |9l . For a subgroup 
T of Q n , we define the subset 

T w = {MeT: wt(M) = w}, 

where wt(M) is the weight of the n-fold Pauli operator M. 
Suppose S is a stabilizer subgroup of the n-qubit Pauli group 
Q n , and Af(S) is its normalize! - group. Let 



A{z) 



w=0 



A w z w , B{z) = Y. B ^ 



w=0 



where 



a w =j2 1 = \ s ™\> B ^= i = iw)u 

Then A(z) is the weight enumerator of the stabilizer S, and 
B(z) is the weight enumerator of the normalizer group M(S). 
Shor and Laflamme proposed the quantum MacWilliams iden- 
tity for general quantum codes [26 1, and Gottesman derived 
the same identity for stabilizer codes [9|. This identity states 
that 



B(z) 



-(l + 3z) n A 



1 



1 + 3z 



or 



1 ™ 

B w =-r^- ^2 P w(w',n)A v 



(3) 
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for w = 0, • • • ,n, where 

w 



u=0 



w \ In — w 
u I \ w — u 



is the Krawtchouk polynomial [19|. This quantum 
MacWilliams identity is exactly the same as the MacWilliams 
identity in ([T]) for a classical linear quaternary code with 
q = 4. 

III. The Dual of an EAQEC Code 

Suppose the sender (Alice) and the receiver (Bob) 
share c pairs of maximally-entangled states = 
^= (|00) + 1 11)). Assume that Alice sends a fc-qubit state \(j>) 
to Bob by using an [[n, k, d; c}] EAQEC code with a Clifford 

is 



encoder U. Then the encoded state \ip} AB 

|^S = (f/ A 0/ BW|^ 0( |^ 



\AB\ 



<£> |0> 



(gm — k— c\ 



where the superscript A or B indicates that the operator acts 
on the qubits of Alice or Bob, respectively. We omit the 
Pauli operators (of each stabilizer generator) on Bob's side, 
and obtained the so-called simplified stabilizer generators: let 
gj = UZjW and hj — UXjW for j = 1, • ■ • , n. Then gi and 
hi are symplectic partners for i = 1, • • • , k + c, meaning that 
they anticommute with each other and they commute with all 
other partners gj and hj, where j ^ i. The simplified stabilizer 
group S' of the encoded state is 

S = (9k+i, ■ ■ ■ ,9k+c, h k+ i, • • • , h k+c , g k + c +i, • • • , g n }- 

Since the generators in S' are simplified, S' is non-abelian. 
The symplectic subgroup of S' is 



Ss — (<7fc+i> " ' > 9k+c, hk+i, 



hk 



+cj 



and the isotropic subgroup of S' is Sj — (g k + c +i, •'• , 9n)- 
Let S' — Ss x Si so that S' is generated by the generators of 
Ss and Si. The logical subgroup C of the encoded state is 

£ = ' ' ' ,9k, hi,-- - ,h k ). 

We begin our exposition of the dual of an EAQEC code 
by considering an entanglement-assisted quantum code with 
logical qubits and ebits alone (the code has no ancillas and 
c = n — k). We call such a code a maximal-entanglement 
EAQEC code. In this case, Si is the trivial group that 
contains only the identity, and the simplified stabilizer group 
is Ss- These EAQEC codes with maximal entanglement are 
in some sense similar to classical codes because they have 
no degeneracy lfT8l . |23'|. Observe that the normalizer of the 
simplified stabilizer Ss is equivalent to the logical group C 
and vice versa: 

M(S S ) = C, N{C) = S s . 

In other words, if the simplified check matrix and logical 
matrix are H and L, respectively, we have 

HA 2n L T = ( 



where 



A 2 „ = 



Or, 



J(2rx2k), 



In 



Oixj is an i x j zero matrix, and I nxn is an n-dimensional 
identity matrix. Thus Ss and C are dual to each other with 
respect to the symplectic product. If we interchange the roles 
of Ss and C, they define another EAQEC code, which is the 
dual of the original EAQEC code. 

Definition 1. The dual of an [[n, k, d;n — k]} EAQEC code, 
defined by a simplified stabilizer group Ss and a logical 
group C, is an [[n, n — k,d';k]] EAQEC code, defined by 
the simplified stabilizer group C and the logical group Ss for 
some minimum distance d! . 

Equivalently, consider that a Clifford unitary encoding op- 
erator acting only on ebits and logical qubits is completely 
specified by the logical operators and the stabilizer operators. 
So to obtain the dual of an EAQEC code, we can just substitute 
all of the information qubits with halves of ebits and vice 
versa. Note that the minimum distance d! of the dual code 
can be determined by the MacWilliams identity for EAQEC 
codes, which will be derived in the next section. 

As an example, consider the class of [[n, 1, n; n — 1]] 
entanglement-assisted repetition codes for n odd [ 18 1. The two 
logical operators for the one logical qubit of this code are as 
follows: 

X X X • • • X X , . , 

(4) 



Z Z Z 



Z Z 



These logical operators are directly imported from the gener- 
ator matrix in The stabilizer generators are as follows: 



Z 


z 


I ■ 


■ I 


I 


I 


z 


z ■ 


■ I 


I 


I 


I 


I ■ 


■ z 


z 


X 


X 


I ■ 


■ I 


I 


I 


X 


X ■ 


■ I 


I 


I 


I 


I ■ 


■ X 


X 



(5) 



One can determine the symplectic pairs by performing a sym- 
plectic Gram-Schmidt orthogonalization of the above operators 
l28l . If we interchange the roles of the stabilizer subgroup and 
the logical operator subgroup, we obtain an [[n, n— 1, 2; 1]] 
accumulator code. 

To make this more precise, consider the encoding circuits 
in Figure [T| The circuit in Figure [TJa) is the encoder of a 
[[9, 1, 9; 8]] entanglement-assisted repetition code. Swapping 
the information qubit for an ebit and all of the ebits for 
information qubits gives the encoder of Figure |TJb), which 
encodes a [[9, 8, 2; 1]] entanglement-assisted accumulator code. 
To illustrate that the circuit is working as expected, let us 
consider it acting on the first five qubits only (just to simplify 
the analysis). Inputting the following two operators at the 
information qubit slot 

X I I I I 
Z I I I I 
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info 
ebit 
ebit 
ebit 
ebit 
ebit 
ebit 
ebit 
ebit 



ebit 
info 
info 
info 
info 
info 
info 
info 
info 



(a) 



(b) 



Fig. 1. (a) The encoder for a [[9, 1,9; 8]] entanglement-assisted repetition code consists of a periodic cascade of CNOT gates. The encoder for arbitrary 
[[n, 1, n; n — 1]] entanglement-assisted repetition codes extends naturally from this design. We can implement these encoders with a simple quantum shift- 
register circuit that uses only one memory qubit 1271 . (b) Considering the circuit in (a) but changing the information qubit to an ebit and all of the ebits 
to information qubits gives the encoder for the dual of the entanglement-assisted repetition code, namely, the [[9, 8, 2; 1]] entanglement-assisted accumulator 
code. This circuit naturally extends to encode the [[n, n — 1,2; 1]] entanglement-assisted accumulator codes. Simple variations of the above circuit encode 
the even n repetition and accumulator EAQEC codes, and we discuss them in Section |VI-B| 



gives the following two logical operators: 



X 

Z 



X 

z 



X 

z 



X 

z 



X 

z 



and these operators match the form of the logical operators for 
the repetition code in Q. Inputting the following operators at 
the ebit slots 



/ 


X 


I 


I 


I 


/ 


z 


I 


I 


I 


/ 


I 


X 


I 


I 


/ 


I 


z 


I 


I 


/ 


I 


I 


X 


I 


/ 


I 


I 


z 


I 


/ 


I 


I 


I 


X 


/ 


I 


I 


I 


z 


operators 






X 


X 


I 


I 


I 


I 


z 


z 


z 


z 


I 


X 


X 


X 


X 


z 


z 


I 


I 


I 


I 


I 


X 


X 


I 


I 


I 


I 


z 


z 


I 


I 


I 


X 


X 


I 


I 


z 


z 


I 



which we can transform by row operations to be the same as 
the operators in Q. 

To extend the notion of a dual code to arbitrary 
entanglement-assisted codes, the approach is to exchange all 
of the information qubits with ebits and vice versa, but do 
nothing with the ancillas. Observe that Af(Ss X Sj) = CxSi 
and J\f{C x Si) — Ss x Si. Thus Ss x Si and C x Si are 
dual to each other. The ancillas have little to do with duality 
in the sense above. That is, the isotropic subgroup Si only 
plays a role in error correcting power. 

Definition 2. The dual of an [[n, k, d; c]] EAQEC code, defined 
by a simplified stabilizer group S' = Ss x Si and a logical 
group C, is the [[n, c, d'; k]] EAQEC code with C x Si being 



the simplified stabilizer group and Ss being the logical group 
for some minimum distance d' . 

We discuss how to determine the minimum distance d! of 
the dual code in the next section. 



IV. The MacWilliams Identity 

We now derive the MacWilliams identity for general 
EAQEC codes. Suppose that the weight enumerators of S' = 
Ss x Si and C x Si are respectively A(z) = X)"=o A w z w 



and B(z) = £™ =0 B w , 



2. 1 = l^|, 

EES', 



w , where 

B v , = 



E 



l = \(£xSi) 



The following lemma is similar to (7.17) and (7.18) in Gottes- 
man's thesis |9l . 



Lemma 3. For an operator E e Q n , 

\E@M _ / 

0, 



E 

MeS' 



if E G CxSi; 
if E 4 C x Si, 



(6) 



where 



EQM 



0, if E and M commute; 

1, if E and M anti-commute. 



Proof: We only give the concept of the proof and a more 
rigorous proof can be provided. If E G C x Si, the result is 
obvious. If E 4_ C x Si, E anti-commutes with half of the 
elements in S 1 . ■ 
Armed with Lemma [3] and Gottesman's proof of the 
MacWilliams identity for stabilizer codes [9|, we have the 
following theorem. 

Theorem 4. If the weight enumerators of S' — Ss x Si and 

C x Si are A(z) and B(z), respectively, then we have the 
MacWilliams identity for EAQEC codes: 



B{z) 



k l+ZzTA {i 



1 - z 

~Tz 
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or 

1 - 

B w =7—77 ^2 p w( w '> n ) A w, for w = 0, ■ • • ,n. (7) 

' ' w'=0 

The Mac Williams identity in (|7]i for EAQEC codes is ex- 
actly the same as the Mac Williams identity in ^ for quantum 
codes. In the case of standard stabilizer codes (c = 0), S$ is 
the trivial group, and Si and C x Si are dual to each other. 
In conclusion, the quantum Mac Williams identity describes a 
relation between the weight enumerators of a group S and its 
normalizer group Af(S). 

In general, the minimum distance d of an EAQEC code 
is the minimum weight of an element in M(Ss x Si)\Si 
iTOl . Suppose C(z) is the weight enumerator of Si. Given 
the weight enumerator B{z) of C x Si and the minimum 
distance d of the EAQEC code, we have C w = B w for w E 
{0, • • • , d — 1}. Equivalently, d is the smallest positive integer 
w such that B w — C w > 0. From the MacWilliams identity, 
we can obtain the weight enumerator A(z) of Ss X 5/. Then 
cf is the smallest positive integer w such that A w — C w > 0. 

In the case of maximal-entanglement EAQEC codes, Si is 
the trivial group and there is no degeneracy. If we exchange 
the roles of Ss and £, we obtain an [[ri, n — k,d'; k]] EAQEC 
code, which is the dual of an [[n, k,d;n — k}) EAQEC code. 
The minimum distance d' of this [[n, n—k, d'; k]] EAQEC code 
is the minimum weight of a nontrivial element in M(C) = Ss- 
Thus d! can be determined from the MacWilliams identity and 
the weight enumerator of the [[n,k,d;n — k]] EAQEC code, 
as in the following example. 

Example 1. The dual of the [[n, 1, n; n—1)] repetition code is 
the [[n, n—1,2; 1]] accumulator code whenever n is odd. The 
coefficients of B(z) for the odd-n [[n, l,n\n — 1]] repetition 
code are B/ n -\ — (1,0, 0,--- ,0,3). Using the MacWilliams 
identity, we obtain the weight enumerators of these dual 
EAQEC codes: 

A {3) =(1,0,9,6), 

B (3 ) =(1,0,0,3), 

A (5) =(1,0,30,60,105,60), 

B (5) =(1,0,0,0,0,3), 

A {7) =(1, 0, 63, 210, 735, 1260, 1281, 546), 

B {7) =(1,0,0,0,0,0,0,3), 

A (9) =(1, 0, 108, 504, 2646, 7560, 15372, 
19656,14769,4920), 

B {9) =(1,0,0,0,0,0,0,0,0,3), 
A (11) =(1, 0, 165, 990, 6930, 27720, 84546, 

180180, 270765, 270600, 162393, 44286), 
B ( u) =(1,0,0,0,0,0,0,0,0,0,0,3) 

We have also verified the correctness of the above weight 
enumerators by counting the logical operators with a computer 
program. 



V. Linear Programming Bounds for EAQEC Codes 

The significance of the MacWilliams identities is that linear 
programming techniques can be applied to find upper bounds 
on the minimum distance of EAQEC codes. We first pro- 
vide the linear programming bound for maximal-entanglement 
quantum codes, which is similar to the case of classical codes, 
since there is no degeneracy in this case. Then we proceed to 
the general case. 

A. The Linear Programming Bounds for Maximal- 
Entanglement Quantum Codes 

The minimum distance of an EAQEC code (or a standard 
stabilizer code) is the minimum weight of any element in (C x 
Si)\Si. Suppose the weight enumerators of Ss x Si and C x 
Si are A(z) and B(z), respectively. For c = or c = n— r, the 
minimum distance can be determined from the coefficients of 
the weight enumerators A(z) and B(z), which are constrained 
by the MacWilliams identities. This gives the constraints we 
need in the integer programming problem. 

For an [[n, k,d;n— k]] EAQEC code, it must have B w = 
for w = 1, • ■ ■ , d — 1. If we cannot find any solutions to an 
integer program with the following constraints: 

A = B = 1; 

A w > 0, B w > 0, for w = 1, • • • , n; 
A w < \S S \,B W < |£|, for w = 1, • • • ,n; 

n n 

^2A W = \S S \,^2B W = \C\; 

1 " 

B w = t— -, } P w (w', n)A w , , for w = 0, • • • , n; 

B w = 0, for w = 1, • • • , d — 1; 

for a certain d, this result implies that there is no [[n, k,d;n — 
k]] EAQEC code. If d* is the smallest of such d's, then d* — 1 is 
an upper bound on the minimum distance of an [[n, k, d; n—k]] 
EAQEC code. This bound is called the linear programming 
bound for EAQEC codes with maximal entanglement. If we 
change the constraint B w = 0, for w = 1, • • • , d — 1 with 

A w =0, for w = 1, • • • , d — 1, 

this gives the linear programming bound on the minimum 
distance of the dual code [[n, n — k,d; k]]. 

Example 2. Consider the [[8, 3, 5; 5]] EAQEC codes from the 
random optimization algorithm in lfl8l . The linear program- 
ming bound shows that there is no [[8, 3, d; 5]] EAQEC code 
with d > 5, and thus the [[8, 3, 5; 5]] code is optimal]^] 

Example 3. Consider the [[15, 7, 6; 8]] EAQEC codes from the 
random optimization algorithm in lfl8l . The linear program- 
ming bound shows that no [[15, 7, d; 8]] EAQEC code with 
d > 7 exists; however, it does not rule out the existence of the 
[[15, 7, 7; 8]] code. 

3 We used the optimization software LINGO from LINDO Systems to solve 
the integer programming problem for these examples. 
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B. The Linear Programming Bound for Non-Maximal- 
Entanglement Quantum Codes 

Now we consider the case < c < n—k, and Si, Ss are not 
trivial. Since the minimum distance of an EAQEC code is the 
minimum weight of any element in (£ x Si)\Si, constraints 
on the weight enumerator of 5/ are required. Observe that 

JV{Si) =LxS s xSi. 

If we view the group (C x Ss) as a logical group of a 
standard stabilizer code with a stabilizer group Si, this is 
just the dual relation of the standard stabilizer codes. Let the 
weight enumerator of C x Ss X Si be D(z). Then C(z) and 
D(z) are related by the Mac Williams identity for standard 
stabilizer codes. The minimum distance of the EAQEC code 
is the minimum nonzero integer w such that B w — C w > 0. 
With constraints on B m 's and C w '%, we can find the linear 
programming bound on the minimum distance of the EAQEC 
code. To sum up, we have the following constraints: 

A) = B Q = C = D Q = 1; 

A w > 0, B w > 0, C w > 0, D w > 0, for w = 1, • ■ ■ , n; 
A w < \S S x Sj\,B w <\Cx Si\, for w = l,--- ,n; 
C w < \Si\,D w < \C x Ss x 5/|, for w — 1, • • • ,n; 
D w > A W ,D W > B W ,D W > C w , for w = 1, • • • , n; 
A w > C W ,B W > C w , for w = 1, • • • , n; 

n n 

Au, = 1 5 S x 5/|, B ™ = \ c x 5 'l> 

iu= 

n n 

J2 C ™ = \ S i\> E A " = l>C x 5 S x 5/|; 



1 - 

S w = — — - } P w (w', n)A w >, for w = 0, • • ■ ,n; 

I^x5/| ^ Q 

1 ™ 

^ = jp-r ^ P w (w',ri)C w >, for w = 0, • • • , n; 
= C^,, for w = 1, • • ■ , d— 1; 

If we cannot find any solutions to the integer programming 
problem with the above constraints for a certain d, this result 
implies that there is no [[n, k, d; c]] EAQEC code. If d* is the 
smallest of such gTs, then d* — 1 is an upper bound on the 
minimum distance of an [[n, k, d; c]] EAQEC code. This is the 
linear programming bound for EAQEC codes with < c < 
n — fc. If we change the constraints 

B w = C w , for w = 1, • • • , d — 1, 

with 

A w — C w , for w = 1, • • • , d — 1, 

this gives the constraints of the linear programming bound on 
the minimum distance of the [[n, c, d; k]] dual code. 

Example 4. The linear programming bound of the [[7, 1;4]] 
code is d < 6, which is the same as the singleton bound. Hence 
the [[7, 1,5; 4]] EAQEC code from the encoding optimization 



algorithm in Example 4 of lfl~8l does not achieve the linear 
programming bound. It follows that the [[7, 1, 5; 5]] EAQEC 
code does not achieve the upper bound either. 

Example 5. The linear programming bound of the [[8, 3; 3]] 
EAQEC code is d < 4, which improves the singleton bound: 
d < 5. Hence the [[8, 3, 4; 3]] EAQEC code from the encoding 
optimization algorithm in Example 8 of lfl~8l is optimal. On 
the other hand, the linear programming bound of the [[8, 3; 4]] 
EAQEC code is d < 5, which is the same as the singleton 
bound. Hence the [[8, 3, 4; 4]] EAQEC code does not achieve 
the upper bound. 

Example 6. The linear programming bounds of [[9, l;c]] 
EAQEC codes for c = 3, 4, 5 are d < 8; however, the singleton 
bounds are d < 6, d < 7, and d < 7, respectively. 

From these three examples, we can determine that the 
linear programming bound might be better or worse than the 
singleton bound when < c < n — k. 

VI. Bounds on EAQEC Codes 

This section establishes a table of upper and lower bounds 
on the minimum distance of maximal-entanglement EAQEC 
codes for n < 15. We begin by discussing the existence of 
arbitrary EAQEC codes, followed by some specific EAQEC 
code constructions. 

The existence of an [[n, k, d\] stabilizer code implies the 
existence of an [[n, k, d! > d; c]] EAQEC code, since we 
can replace ancilla qubits with ebits and then optimize the 
encoding operator [18]. Therefore, the lower bound on the 
minimum distance of standard codes [8| can be applied here. 
Similarly, the existence of an [[n, k, d; c]] EAQEC code where 
c < n — k implies the existence of an [[n,k,d' > d; c' > c]] 
EAQEC code. This basically establishes the existence of many 
[[ti, fc; c]] EAQEC codes. 

A. Gilbert-Varshamov Bound for EAQEC Codes 

Consider the stabilizer group 5 of an [[n, k, d; c]) EAQEC 
code, which is a subgroup of the Pauli group Q n + C . We 
consider only the error operators in the group Q n because 
the entanglement-assisted paradigm assumes that the ebits on 
Bob's side of the channel are not subject to errors. Now the 
EAQEC code is defined in an (ti + c)-qubit space, but only 
the first n qubits suffer from errors. Following the argument 
of the quantum Gilbert-Varshamov bound [4), we can get the 
Gilbert-Varshamov bound for EAQEC codes. However, we 
will show that there are maximal-entanglement EAQEC codes 
with minimum distance higher than the value predicted by this 
bound for n < 15. 

Theorem 5. Given n, d, c, let 

r / 2™+ c \ 

The Gilbert-Varshamov bound for EAQEC codes states that if 
< k < n — c, then there exists an [[n, fc, d; c\] EAQEC code. 
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Equivalently, 



For example, when n = 6, 



E3 J 

3=0 



2 fc > 2 



n+c 



B. Maximal-Entanglement EAQEC Repetition and Accumula- 
tor Codes for Even n 

Lai and Brun proposed a construction of [[n, l,n;n — 1]] 
entanglement-assisted repetition codes for n odd in ifTHl . In the 
case of even n, that construction gives a series of [[n, 0,n;n — 
2]] EA repetition codes with no information qubits. 

In this section, we construct [[n, 1, n — l;n — 1]] EA 
repetition codes for n even. We obtained these codes both 
with the techniques from Ref. |[T8l and by realizing that the 
encoding circuit in Figure [T] can encode them. The duals of 
these codes are the [[n,n — 1,1; 1]] EA accumulator codes 
for even n. Theorems [7] and [8] show that both of these code 
constructions are optimal, in the sense that [[n, 1, n; n— 1]] and 
[[n, n— 1,2; 1]] EAQEC codes do not exist for even n. Thus, 
the results here complete our understanding of the dual classes 
of entanglement-assisted repetition and accumulator codes for 
arbitrary n. 

Theorem 6. There are [[n, 1, n — 1; n — 1]] entanglement- 
assisted repetition codes for n even. The duals of these codes 
are the [[n,n — 1,1; 1]] entanglement-assisted accumulator 
codes. 

Proof: Suppose Hf n _i) is an (n — 2) x (n — 1) parity- 
check matrix of a classical [n — 1, 1, n — 1] binary repetition 
code: 

"1 1 •••0 0' 
1 1 ••• 



H 



(n-l) 







1 1 



We define two (n — 1) x n matrices 
"0 



H, = 



■ff(n-l) 





1 1 



and 



Ho 



H 



(n-l) 



1 1 



1 
0" 



1 1 



Consider a simplified check matrix of the form 

H' = 



O Hi 
H 2 O 



Consider the matrix We have that 

if i = j for j — 1 , • • • , n — 1 , or 



1 



i = j -2 for j = 3,--- ,n-2 
0, else. 



h%ht — 



10 10 
10 10 
10 
1 
1 

Thus the number of symplectic pairs in H' is as follows IT71 : 
1 



rank(iJ'A J ff' / ) = rank( J ffiff 2 i ) 



1. 



The simplified logical 


matrix 


is 




L' = 


' 00- •• 


00 


11-1 " 


11... 


10 


00 • • • 



which implies the minimum distance is n — 1. Therefore, H' 
and V define an [[n, 1, n— 1; n— 1]] EAQEC code. One obtains 
the dual [[n, n — 1, 1; 1]] codes simply by swapping the roles 
of the logical matrix and the simplified check matrix. 

■ 

This completes the family of entanglement-assisted repetition 
and accumulator codes for any n. The encoding circuit of 
Figure [T] encodes these even-n repetition codes with the 
exception that the last qubit is removed, the last CNOT in 
the first string does not act, and the last CNOT in the second 
string does not act. 

Example 7. The coefficients of B(z) of the even-n [[n, 1, n — 
1; n — 1]] repetition code are Br n \ = (1, 0, • • ■ ,0, 1, 2). Using 
the MacWilliams identity, we obtain the weight enumerators 
of these dual even-n EAQEC codes: 



(4) 
(4) 
(6) 



A (s) = 

-8(8) = 

-4(10) = 



B 



(10) 
1(12) 



(1,1,15,27,20), 
(1,0,0,1,2), 

(1,1,40,130,305,365,182), 
(1,0,0,0,0,1,2), 

(1, 1, 77, 357, 1435, 3395, 5103, 4375, 1640), 
(1,0,0,0,0,0,0,1,2), 
(1, 1, 126, 756, 4326, 15246, 38304, 65604, 
73809,49209,14762), 
(1,0,0,0,0,0,0,0,0,1,2), 
(1, 1, 187, 1375, 10230, 47850, 168630, 432894, 
811965, 1082565, 974303, 531443, 132860), 
(1,0,0,0,0,0,0,0,0,0,0,1,2). 



B (12 ) = 

The even-n [[n, l,n — l;n — 1]] repetition codes do not 
saturate the quantum singleton bound or the linear program- 
ming bounds. Were an even-n [[n, 1, n; n— 1]] code to exist, it 
would have a weight enumerator B(z) with Bq = 1, B n = 3, 
and C w — for w ^ 0, n. The weight enumerator of its dual 
would also have the the coefficients 



which are positive integers for w = 0, • • ■ , n. It would only be 
able to correct up to I ^i^J channel qubit errors, which is the 
same number of errors that our even-n repetition codes can 
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correct. We prove below that even-n [[n, 1, n; n — 1]] codes do 
not exist, and thus our even-n repetition codes from Theorem[6] 
are optimal. 

Theorem 7. There is no [[n, 1, n; n — 1]] EAQEC code for n 
even. 

Proof: We prove it by contradiction. Suppose there is an 
[[n, 1, n; n — 1]] EAQEC code for n even with a 2 x 2n logical 
matrix 







" u 1 


v 1 ' 






. u2 


v 2 _ 



where u 1 , u 2 , v 1 , and v 2 are binary vectors of length n. These 
vectors should satisfy the following condition in order for the 
above matrix to be a valid logical matrix: 

u 1 ■ v 2 + u 2 ■ v 1 = 1 mod 2. 

Let gw( ) be the "general weight" function defined by 

gw(u\v) = !> 

i:Ui — 1 or Vi=l 

where ui denotes the i th bit of the binary n-tuple u. The above 
binary vectors should satisfy the further constraints 

gw (u 1 1 v 1 ) = gw (u 2 \v 2 ) = gw (u 1 + u 2 {v 1 + v 2 ) = n, 

in order for the code to have distance n as claimed. 

We now exploit the above constraints in order to obtain a 
contradiction. We first partition the first row of the matrix into 
subsets A, B, and C of X, Y, and Z operators, respectively. 
There should not be any identity operators in the first row in 
order for the code to have distance n. Up to permutations on 
the qubits (under which the distance is invariant), the logical 
matrix has the following form: 











where 1 is a vector of all ones, is a vector of all zeros, and 
we have split up the vector (u 2 \v 2 ^ into different components 
corresponding to the subsets A, B, and C. Consider the vector 



Suppose that a component 



0. Then 



should 



equal 1 so that the code's distance is not less than n. Now 
suppose that (u\) . = 1. Then (v\) . should also equal 1 so 
that the code's distance is not less than n. Otherwise, we could 
add ( 1 1 0) to ((u\) . | (Va) ■) and obtain (0|0) as a component 
of another logical operator, and such a result would imply that 
the code's distance is less than n. These steps imply that the 
logical matrix should have the following form: 



1 

,,2 



1 

,,2 







•c 





1 



1 

,.2 



"C 



Similar 
equal 1, 
form: 



reasoning with Vq and Uq implies that 



should 



and the logical matrix should then have the following 



1 

,,2 



1 

,,2 





1 



1 

,.2 



1 

•4 



Finally, consider the vector u 2 B . Suppose that a component 
(Ug) . = 1. Then (v%) i should equal so that the code's 
distance is not less than n. Otherwise, we could add (1|1) to 



(( M s)i I anc ^ obtain (0 1 0) as a component of another 

logical operator, and such a result would imply that the code's 
distance is less than n. Now suppose that . = 0. Then 
(Vg) . should equal 1, by reasoning similar to the above. Thus, 
the logical matrix should have the following form in order for 
the code's distance to be equal to n: 



1 





1 



1 



1 



(8) 



where u 2 B is the binary complement of u 2 B . Now, the symplec- 
tic product of the above two vectors is 

{\A\ + \B\ + |C|)mod2 = nmod2 = 0, 

which contradicts the assumption that the original matrix is a 
valid logical matrix. ■ 

Theorem 8. There is no [[n, n — 1, 2; 1]] EAQEC code for n 
even. 



Proof: We prove the theorem by contradiction, in a 
fashion similar to the proof of the previous theorem. Suppose 
there is an [[n,n — 1, 2; 1]] EAQEC code for n even, and 
suppose its 2 x 2n simplified check matrix [Hx\Hz] has the 
form 









v 1 ' 






. « 2 


v 2 _ 



where u , u 2 , v , and v 2 are binary vectors of length n. These 
vectors should satisfy the following condition in order for the 
above matrix to be a simplified check matrix of a maximal- 
entanglement EAQEC code with one ebit: 



u ■ v -\- u 



1 mod 2. 



We now partition the first row of the simplified check ma- 
trix into subsets A, B, and C of X, Y, and Z operators, 
respectively. Up to permutations on the qubits (under which 
the distance is invariant), the simplified check matrix has the 
following form: 



1 



1 





u 2 c 





i< 2 



1 



1 



where we have split up the vector (u 2 \v 2 ) into different 
components corresponding to the subsets A, B, and C. The 
code has a minimum distance two by assumption and is 
non-degenerate because it is a maximal-entanglement EAQEC 
code. Therefore, no column of the above matrix should be 
equal to the all-zeros vector. Were it not so, then the code 
would not be able to detect every single-qubit X or Z error 
and would not have distance two as claimed. These constraints 
restrict the simplified check matrix to have the following form: 



1 

4 



i 





1 



1 

"I 



1 



Also, no column of the entrywise addition of the matrix to 
the left of the vertical bar and the matrix to the right should 
be equal to the all-zeros vector. Were it not so, then the code 
would not be able to detect every single-qubit Y error and 
would not have distance two as claimed. These constraints 
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further restrict the simplified check matrix to be as follows: 



1 

< 2 y 



1 

,2 





1 



1 

ti% 



(9) 



u A «• 

Now, the symplectic product of the above two vectors is 

{\A\ + \B\ + |C*|)mod2 = nmod2 = 0, 

which contradicts the assumption that the original matrix is a 
simplified check matrix for a maximal-entanglement EAQEC 
code with one ebit. ■ 
Interestingly, observe that the non-existent logical matrix in 
(|8]l has the same form as the non-existent simplified check 
matrix in |9]). Were either type of code to exist, we would 
expect them to be duals of each other, but they both fail to 
exist because they cannot satisfy the dual constraints imposed 
on them! 

C. Existence of Other EAQEC Codes 

The following theorem is similar to Theorem 6 in Ref. (SJ. 
It shows how to obtain new EAQEC codes from existing ones. 
These results are helpful in our search for lower bounds on the 
minimum distance of maximal-entanglement EAQEC codes. 

Theorem 9. Suppose an [[n, k, d; c]] code exists. Then 

1) An [[n + 1, k, d; c + 1]] code exists. 

2) An [[n, k — 1, d! > d; c + 1]] code exists. 

Proof: 

1) Suppose H = [Hx\Hz] is a simplified check matrix of 
an [[n, k, d; c]] code. Then the simplified check matrix 



H' 



defines an [[n + 1, k, d: c+1]] code. We have the stabilizer 
group (S<S)I)U{X n+ i,Z n+ i}, where (5® J) = {E®I : 
EeS}. 

2) It is obtained by moving a symplectic pair from the 
logical group to the stabilizer group. 



D. The Plotkin bound for EAQEC Codes 

The Plotkin bound for EAQEC codes is similar to the 
Plotkin bound for classical codes |[T9ll . It is again helpful 
in our efforts to bound the minimum distance of maximal- 
entanglement EAQEC codes. 

Theorem 10. The Plotkin bound for any [[n, k, d; c]] code is 

3nM 



" 00 • • • 








00- •• 





1 " 


00- •• 





1 


00- •■ 






















H x 























_ 



d< 



4(M - 1) 



where M 



y21, 



Proof: The proof is based on the proof of the classical 
Plotkin bound in ifljO- Let M = 2 2k be the number of logical 
operators in the code. We bound the quantity J2 U ug£ wt(u-v) 



in two different ways. First, we lower bound it. There are M 
choices for u, and for each choice of u, there are M — 1 
choices for v. Furthermore, for a code of minimum distance 
d, wt(u ■ v) > d for any u ^ v. So the following lower bound 
holds 

M(M-l)d< E wt(u-u) 

— 2_, Wt ( U ' v )- 

The equality holds because wt(u • v) = if u = v. Now we 
obtain an upper bound on the quantity. We form an M x n 
matrix whose rows are the elements in the logical group C. 
Let m\, m 3 2 , m 3 3 , and m\ be the number of /, X, Y, and 
Z operators in column j of this matrix, respectively. So the 
equality 2~2t=i m \ = M holds for all j e {1, • • • , n}. Each 
choice of a particular Pauli operator and some other Pauli 
operator in the same column contributes exactly 2 to the sum 
vec wt ( M ' v )- Thus, the first equality below holds for this 



reason, and the second holds by applying 2~2i=i '■ 

n 4 

u,v££ j—1 l—l 



M: 



= t ( M2 


-EK') 2 




1=1 


< ± (m* 








= — M . 
4 







The first inequality follows by applying 2~2t=i TO / /4 = M/4 
and convexity of the squaring function: 

/ 4 \ 2 4 2 

(m/4) 2 = E m '/ 4 ^E( m / 4 - 



1=1 



Combining the lower and upper bounds gives us the Plotkin 
bound. ■ 
Since the proof is independent of the number of ebits c, the 
Plotkin bound applies to arbitrary EAQEC codes. However, 
note that c does not appear in the bound, and consequently, 
this bound best describes the characteristics of maximal- 
entanglement EAQEC codes. However, for large k, the bound 
is approximately |n. Hence, this bound is useful only for 
small values of k. 

E. Table of Lower and Upper Bounds on the Minimum Dis- 
tance of Maximal-Entanglement EAQEC Codes 

Combining all the results in this and the previous section, 
we establish Table [j] which details lower and upper bounds 
on the minimum distance of maximal-entanglement EAQEC 
codes with length n < 15. The lower bounds in Table [I] are 
slightly higher than the Gilbert- Varshamov bounds for n < 15. 
The matching upper and lower bounds for k — 1 are from the 
family of entanglement-assisted repetition codes. The upper 
bounds for n < 15 and k > 2 are from the linear programming 
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in 



7; 



bound, which is generally tighter than the singleton bound: 

n - k + c > 2(d - 1), 
and the Hamming bound for non-degenerate EAQEC codes 

m-. 

J^Y") < 2 "- fc+c . 
j'=o 

The Plotkin bound and the linear programming bound match 
for k < 2 and n < 15. For k = 3 and n = 
4,5,6,9,10,11,13,14,15, they also match. For k > 3, the 
Plotkin bound is not as tight as the linear programming bound, 
the singleton bound, or the Hamming bound. 

The following codes are from |fl8l : [[7, 2, 5; 5]] 
Example 4; [[8, 3, 5; 5]], [[13, 3, 9; 10]] in Example 
[[15, 7, 6; 8]], [[15, 9, 5; 6]] [[15, 8, 6; 7]] in Example 
[[7, 3, 4; 4]], [[8, 2, 6; 6]], [[10, 3, 6; 7]] in Table 10. 

The following codes are from fl29]|: [[6, 2, 4; 5]], [[8, 2, 4 
[[9, 6, 3; 3]], [[10, 6, 4; 4]], [[10, 7, 3; 3]], [[11,8,3: 
[[12, 6, 5; 6]], [[12, 7, 4; 5]], [[12, 9, 3; 3]], [[13, 5, 6; 8]], 
[[13, 9, 4; 4]], [[13, 10, 3; 3]], [[14, 11, 3; 3]], [[15, 4, 8; 11]], 
[[15, 10, 4; 5]]. 

From part 1) of Theorem [9] we have [[11, 6, 4; 5]], 
[[12, 8, 3; 4]], [[14, 3, 9; 11]], [[14, 9, 4; 5]], [[15, 11, 3; 4]] codes. 
From part 2) of Theorem |9j we have [[10, 5, 4; 6]], 
[[13, 8, 4; 5]], [[14, 10, 3; 4]] codes. From part 2) of Theorem 
[9] and odd-n EA accumulator codes, we have the codes of 
minimum distance 2. From Theorem [8] and the even-n EA 
accumulator codes, we have the codes of minimum distance 
1. 

We used MAGMA ||30) to find the optimal [[n, k, d]] codes, 
and then applied the optimization algorithm in Ref. ifTSI to 
obtain the other lower bounds. 



;4]], 
;3]], 



VII. The Weight Enumerator Bound on the Block 
Error Probability under Maximum Likelihood 
Decoding 

Since maximal-entanglement codes bear many similarities 
to classical codes, the block error probability when trans- 
mitting coded quantum information through the depolarizing 
channel can be upper bounded with the weight enumerator of 
a particular maximal-entanglement EAQEC code (similarly to 
the case for classical codes \24\, [25 1). This "weight enumer- 
ator bound" gives an idea of the performance of maximum- 
likelihood decoding of an arbitrary maximal-entanglement 
EAQEC code. We can also determine the expected perfor- 
mance when decoding a random EAQEC code with a max- 
imum likelihood decoding rule. Below, we determine these 
bounds and plot them for the maximal-entanglement repetition 
and accumulator EAQEC codes. The result is that these codes 
perform comparably to a random EA code with respect to this 
upper bound. 

Theorem 11. Suppose that the sender transmits an 
[[n, k;n — k)] maximal-entanglement EAQEC code over a 
depolarizing channel with parameter p, and furthermore, that 
the receiver decodes this code according to a maximum a 
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TABLE I 

Upper and lower bounds on the minimum distance of any 

[[n, k,d;n - k]] MAXIMAL-ENTANGLEMENT EAQEC CODES. 



posteriori (MAP) decoding rule. Then we have the following 
upper bound on the block error probability Pb'- 



P B < J B( 7 )-1, 



(10) 



where B (z) is the weight enumerator of the maximal- 
entanglement EAQEC code and 7 is the "Bhattacharyya pa- 
rameter" for the depolarizing channel: 

7-2^f(l-p) + |p. 

Proof: Let U be a Clifford encoder for the [[n, k; n — k]] 
maximal-entanglement EAQEC code. The encoded state 
\^) AB is \i>) AB = (U A ® I B ) ® . 
Then Alice transmits her qubits (entangled with Bob's qubits) 
through n independent uses of a depolarizing channel £ where 

£ (p) = (1 - p) p + I (XpX + YpY + ZpZ) , 

and p is the density operator of a single qubit. We assume that 
p < 3/4 because the channel is completely depolarizing when 
p = 3/4. Suppose that an error operator E £ Q n occurs after 
the depolarizing channel, and that s x ,s z are the binary vector 
representations of the error syndrome. Both s x and s z are of 
length (n — k), and Bob observes them by first decoding the 
qubits with a decoding unitary W and then performing Bell 
measurements on the ebits. This implies that 



{U^E) A ®I B \^) AB 

=Z O |0)® ((x sX z : 



®I B )(\<1> + ) AB ) 
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and E = U{Lq®X s Z s for some logical error Lq € Co, 
where Cq is the set of unencoded logical operators. Poulin 
et al. devised a maximum a posteriori decoder for standard 
stabilizer codes [31 iQand we can modify their decoder to be 
a maximum a posteriori decoder for maximal-entanglement 
EAQEC codes Lmap (s x ,s z ), where 

£map(s x ,s z ) = argmaxPr{L|s x ,s z } . 

Lec 

This decoder selects the most likely recovery operation on the 
decoded qubits, given the syndrome information s x and s z . We 
can calculate the above conditional distribution by applying the 
Bayes rule to the joint distribution Pr {L, s x , s z }: 

Pr{L,s x ,s z } 



Pr{L|s x ,s z } 



Y, v Pr{i',s x ,s z }' 



where 



Pr{i,s x ,s z } = Pr{^}| B=a(i8X8 x zs , )£/t 



(i-pY 



(i-pY 



-wt(J3) 



(!) 



Wt(B) 



3(1 -p) 



E=U (L®X* x Z sZ )[/t 
wt(B) 



The distribution Pi" s x , s z } is fixed over all choices 
of L satisfying the constraint P = U (L ® X s " Z s ") W 
where U is the encoding unitary. Since p < 3/4 <^=^ 
p/(3(l—p)) < 1, the best choice of L for the maximum 
a posteriori decoder Lmap(s x ,s z ) is the one that selects a 
recovery operator L _1 = L such that E = U{L®X S ^ Z S *)W 
with the minimum weight. This minimum distance decoder is 
similar to a classical minimum distance decoder. 

Let C be the set of encoded logical operators and S' be the 
set of simplified stabilizer generators. Given Eq £ Co, let 

Q(E ,L ) = {s x ,s z :Pr{Z £ ,s x ,s z } > Pr{L , s x , s z }}. 

We can now bound the probability Pb(Lq) of a block error 
given that the error operator L occurs under this decoding 
scheme: 

P B (L ) =Pr{MAP decoder fails| L occurs} 
= Pr{L MA p (s x ,s z ) ^L } 
= Pr [L -L MAP (s x ,s z )^/} 

= Pr{Zo-L M AP (s x ,s z )eA,V) 
< V Pr{Z -L M Ap(s x ,s z ) =E ] 



Eo&CoV 

^ E E 

E eC \I S *, S *£Q(E ,L ) 



Pr{Z £ ,s x ,s z } 



Since \ P * {L ( f°ff, } > 1 for s x ,s z e Q(E ,L ), we can 



4 Poulin et al. described their decoder as a "maximum-likelihood" decoder 
|31|, but a careful study of it reveals that their decoder should more properly 
be called a maximum a posteriori decoder. 



multiply each term in sum by this factor and then 

Pb(L ) 



< E E JPr{io,s x ,s z }Pr{L £ ,s x ,s z } 

E eC a \I s x ,s^eQ(E ,L ) 



- E E JPr{io,s x ,s z }Pr{Z £o,s x ,s z } 



= E E Jp^{lm}pt{lem} 

E£C\IMGS' 

< E E Jpt{Lm}Pt{lEm}, 

E£C\I M£Q n 

where L = UL Q W, E = UE W, and M 
U(X sX Z s *)W € S'. Observe that 



Meg" 



Jpt{lm\Pt{lEm\ 



MeG" 



Jpt{M}Pt(lELm\ 



M£G n i=l 



= 11 E JptKM^Pt^LUEULUM)^. 
It holds that (L) i (E) i (L) l ^ I if (E)i ^ I and so 



(M) i ee 



£ JPr{(M)i}Pr W<} 



= 2j|(l-p) + |p = 7 . 



Otherwise, 



= 1. 



J] A/Pr { (M)j } Pr |(i)i(E)i(L)j(M)j| 
(M)ieS 
Consequently, 

Pb(L )< £ ^ (iEL) 

E<EC\I 

= E iWt{E) 

Eec\i 
= 5(7) ~ 1. 

Therefore, the probability Pb of a block error is bounded by 
.6(7) — 1 when taking the expectation over all Lq. ■ 
The above theorem is similar to Theorem 7.5 in Ref. Il25l . 
which determines an upper bound on the block error proba- 
bility when transmitting a classical linear code over a binary 
symmetric channel. 

Theorem 12. Suppose that the sender transmits a random 
[[n, k;n — ft]] maximal-entanglement EAQEC code over a 
depolarizing channel with parameter p and furthermore that the 
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receiver decodes this code according to a maximum likelihood 
decoding rule. Let U be the Clifford encoder for this code. 
Then we have the following upper bound on the expected 
block error probability Pb- 

P b =Eu{Pb}< ((1 + 37)" - 1) , (ID 

where 7 is defined in the previous theorem and the expectation 
is with respect to the choice of random code. In particular, if 
the rate k/n satisfies the following upper bound: 

-<l-ilog a (l + 37), 
n 2 

then the error probability decreases exponentially to zero in 
the asymptotic limit. 

Proof: We first establish a method for choosing a random 
maximal-entanglement EAQEC code. A natural method for 
doing so is first to fix a basis of Pauli operators X\, Z\, X%, 
Zi, ■ ■ ■ , X n , Z n , where the first n—k anticommuting pairs 
correspond to the stabilizer operators for the n—k ebits and the 
next k anticommuting pairs correspond to the logical operators 
for the k information qubits. We then select a Clifford unitary 
uniformly at random from the Clifford group (see Section VI- 
A-2 of Ref. IT321 for a relatively straightforward algorithm for 
doing so) and apply it to the above fixed basis. This procedure 
produces 2n encoded operators Xi, Zi, X 2 , Z 2 , ■ ■ ■ , X n , Z n 
that specify a random maximal-entanglement EAQEC code. 

We now need to determine the expected weight enumer- 
ator Ejj {B(z)} = Y^iw=o ^ u {Bw} z w for such a random 
maximal-entanglement code because it would allow us to 



apply Theorem 11 in order to get an upper bound on the 
expected block error probability. Each coefficient Ejj {B w } 
corresponds to the expected number of Pauli operators of 
weight w that belong to the logical operator group of a 
random EA code. Equivalently, it corresponds to the expected 
number of Pauli operators of weight w that commute with 
the entanglement subgroup of a random code. First, let us 
consider Ejj {-Bo}- The identity operator is the only Pauli 
operator with weight zero. It commutes with all operators with 
unit probability. Thus, Ejj {Bo} = 1. Now, let us consider 
Ejj {B w } with w > 1. We first determine the probability that 
a Pauli operator g with non-zero weight commutes with the 
2 (n — k) encoded operators X\, Z\, X 2 , Z 2 , . . • , X n -k, 
Z n _k for a random EA code. To simplify the calculation, 
observe that applying a uniformly random Clifford unitary to 
the operators X\, Z\, X 2 , Z 2 , ■ ■ ■ , X n -k, Zn-k an d then 
determining the probability that a fixed operator g commutes 
with all of them is actually the same as keeping the basis fixed 
and applying a random Clifford to the operator g itself. This 
holds because 



CfC^g±gCfC^ = 



ftfgC±ClgCf = 0. 



Then a uniform distribution on the Clifford unitaries takes this 
operator g to an arbitrary Pauli operator g' , and the distribution 
induced is just the uniform distribution on all of the 2 2n — In- 
qubit Pauli operators not equal to the identity (this reasoning 
is the same as that in Section VI- A- 1 of Ref. [32j). At this 
point, the argument becomes purely combinatorial, and the 



only operators that commute with the above fixed basis are the 
ones with identity acting on the first n — k qubits. Thus, there 
are 2 2k — 1 Pauli operators besides the identity that commute 
with the fixed basis, and we conclude that the probability that 
a fixed Pauli operator g commutes with the random set X\, 
Zi, X 2 , Z 2 , !„_);, Z n -k is 



■,2k 



1 



2 2 ™ - 1 ' 

Now we can calculate the expected number of operators that 
are in the logical subgroup. The number of Pauli operators 
with weight w is (™)3 lu . Consequently, we have 

2 2k 2 

Eu{B w } = 



2 2n _ 1 



which implies 

Eu{B(z)}-Eu{B } 



w—1 
2 2fe _ J 

22^ 
2 2fc 



E v {B w }z v 



3*V 



((1 + 3z) n - 1) ■ 



2 2n _ 1 

Therefore, by exploiting the result in Theorem [TTJ an upper 
bound on the expected block error probability for general 
EAQEC codes with maximal entanglement is 



Eu{Pb} < B(7) - Bo, 
2 2k -1 



2 2n _ I 



((l + 3 7 ) n -l). 



We can drive the expected error probability to be arbitrarily 
low in the large n and k limit by ensuring that 



k < 1 - \ log 2 (1 + 3 7 ) ■ 
n 2 



(12) 



This bound is not as tight as the entanglement-assisted hashing 
bound (the optimal limit), and Figure|2]displays how these two 
bounds differ. ■ 
We can plot the error probability bound in ( [T0| as a function 
of p for specific codes such as the repetition codes or the 
accumulator codes and then compare the results with the 
average error probability bound for a random code. Figure 
[3] provides such plots and compares their performance with a 
random EA code, with respect to these bounds. 

VIII. Discussion 

In this paper, we studied the properties of EAQEC codes. 
We defined the dual code of an EAQEC code and derived 
the Mac Williams identities for EAQEC codes, which are the 
same as for the classical case. Based on these identities, 
we found a linear programming bound on the minimum 
distance of an EAQEC code. We also derived the Plotkin 
bound and the Gilbert- Varshamov bound for EAQEC codes, 
together with several theorems examining the existence of 
EAQEC codes. Combining these results, we provided a table 
of upper and lower bounds on the minimum distance of 
maximal-entanglement EAQEC codes for n < 15. Finally, 
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Maximal-Entanglement Repetition Codes 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Depolarizing rate 



(a) 



Maximal-Entanglement Accumulator Codes 




10" 3 10" 2 10"' 10° 



Depolarizing rate 

(c) 



Random Maximal-Entanglement Code with Rate 1/n 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Depolarizing rate 

(b) 



Random Maximal-Entanglement Code with Rate (n-1)/n 




Depolarizing rate 



(d) 



Fig. 3. The figures plot the weight enumerator bound in (TTJ as a function of the depolarizing parameter p for various finite-length codes, (a) The weight 
enumerator bound for maximal-entanglement repetition codes of length 3 to 12. (b) The expected weight enumerator bound for random rate 1/n maximal- 
entanglement codes of length 3 to 12. (c) The weight enumerator bound for maximal-entanglement accumulator codes of length 3 to 12. (d) The expected weight 
enumerator bound for random rate (n — 1) /n maximal-entanglement codes of length 3 to 12. Observe that the performance of the maximal-entanglement 
repetition and accumulator codes with respect to this upper bound is comparable to the expected performance of random maximal-entanglement codes. 



we determined "weight enumerator bounds" on the block 
error probability when decoding according to a maximum- 
likelihood decoding rule, and we found that the performance 
of maximal-entanglement repetition and accumulator codes is 
comparable to the expected performance of random codes, 
with respect to this upper bound. 

The table of upper and lower bounds on the minimum 
distance of any [[n, k, d}] standard stabilizer codes (c = 0) 
is given in |8|. Similar tables for EAQEC codes with < c < 
n — k can be constructed. 

We proposed a construction of [[n, l,n — l;n — 1]] 
entanglement-assisted repetition codes for n even, which com- 
pletes the family of entanglement-assisted repetition codes for 
any n. These entanglement-assisted repetition codes are the 
optimal codes that encode a single information qubit. We also 
constructed an explicit encoding circuit for these codes. 



Most lower bounds in Table [I] are from the optimization 
algorithm [18|. However, these codes may not be close to 
optimal when n + k becomes large, for the complexity of 
the encoding optimization algorithm increases exponentially 
with n + k, making full optimization impossible. To make 
the bounds in Table [j] tighter, we need to consider other code 
constructions to raise the lower bounds. One possibility is to 
extend the encoding optimization algorithm from the form in 
|fT8l , by starting the optimization from a good EAQEC code 
and adding additional entanglement. From a given [[n, k, d; c]] 
code, we can change any subsets of the n — k — c ancilla 
qubits with ebits and then apply a similar encoding optimiza- 
tion algorithm. From a given "good" [[n, k, d±; c]] code, we 
obtain an optimized [[n,k,d';n — k]} code much faster than 
using the original optimization over an [[n, k, ^2]] code. Since 
we begin with a "good" [[n, k, d\; c]] code, the optimized 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
Depolarizing parameter 



Fig. 2. The figure plots both the entanglement-assisted hashing bound 
1 — 1/2 [H2 (p) + p log 2 31 from Ref . 1221 and the "asymptotic weight 
enumerator bound" from l |12( as a function of the depolarizing parameter. 
The two bounds become close for high depolarizing noise. Interestingly, the 
thresholds of the maximal-entanglement EA turbo codes from Ref. | 23J are 
just shy of the asymptotic weight enumerator bound (see Figures 6(b) and 
7(b) of that paper). 



[[n, k, d'; n— k]] code is expected to approach the upper bound. 
Such a "good" code can be obtained from the construction in 
ifPJl , using optimal classical quaternary codes. 

On the other hand, we proved the non-existence of 
[[n, l,n;n — 1]] or [[n,n — 1, 2; 1]] codes for n even, which 
decreases the upper bound predicted by the linear program- 
ming bound for k = 1 and n even. We plan to explore the 
existence of other [[n, k, d;n — c\] codes to decrease the upper 
bound. 

Shadow enumerators, designed for self-dual codes, give 
additional constraints on the weight enumerator of a stabilizer 
code. However, there is no general form in the case of EAQEC 
codes. 

Finally, consider the possibility of a "self-dual code" 
[[n, n/2, d; n/2}] for n even, such that the dual code is also 
an [[n, n/2, d; n/2]] code with the same weight enumerators. 
That is, B(z) — A(z). We conjecture that such self-dual codes 
exist. If so, the two groups Ss and C may be equivalent up 
to a permutation on the qubits. 
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